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ABSTRACT 

We study the cosmological evolution of scalar fields with arbitrary potentials in the presence of a barotropic 
fluid (matter or radiation) without making any assumption on which term dominates. We determine what 
kind of potentials ^(0) permits a quintessence interpretation of the scalar field <j) and to obtain interesting 
cosmological results. We show that all model dependence is given in terms of A = —V /V only and we study 
all possible asymptotic limits: A approaching zero, a finite constant or infinity. We determine the equation 
of state dynamically for each case. For the first class of potentials, the scalar field quickly dominates the 
universe behaviour, with an inflationary equation of state allowing for a quintessence interpretation. The 
second case gives the extensively studied exponential potential. While in the last case, when A approaches 
infinity, if it does not oscillate then the energy density redshifts faster than the barotropic fluid but if A 
oscillates then the energy density redshift depends on the specific potential. 
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Models with a cosmological constant term, intended as a constant vacuum contribution or as a 
slowly decaying scalar field, have recently received considerable attention for several reasons, both 
theoretical and observational. 

From the theoretical point of view, we have to face the possible conflict between the age of the 
universe in the standard Einstein-de Sitter model and the age of the oldest stars, in globular 
clusters. Estimates of the Hubble expansion parameter from a variety of methods seem to point to 
Hq « 70ibl0 km/s/Mpc (a recent review can be found in [|^; see e.g. [|| for specific projects), leading 
to an expansion age of t;/ ~ 9 it 1 Gyr for a spatially flat universe with null cosmological constant. 
On another hand, the age of globular clusters have been estimated in the range ~ 13 — 15 Gyr @; 
although revised determinations based on the Hypparcos distance scale are lower by approximately 
2 Gyr 

The requirements of structure formation models also suggest a cosmological constant term. Simu- 
lations of structure formation profit from the presence of matter that resists gravitational collapse 
and A CDM models provide a better fit to the observed power spectrum of galaxy clustering then 
does the standard CDM model ||, Q. 

On the observational side, we find direct evidence in recent works on spectral and photometric 
observations on type la supernova that favour eternally expanding models with positive cos- 
mological constant. Statistical fits to several independent astrophysical constraints support these 
results Q. See however for a different explanation to these observations. More indirect evidence 
comes from the observational support for a low matter density universe from X-ray mass estima- 
tions in clusters [^], [11|. In these works, if the nucleosynthesis limits on the baryonic mass are 



to be respected, the total matter that clusters gravitationally is limited to J; 0.3. In such a case, 
a cosmological term would reconcile the low dynamical estimates of the mean mass density with 
total critical density suggested by inflation and the flatness problem. 

Many models with a scalar field playing the role of a decaying cosmological constant have been 
proposed up to now. Some of them are specific models motivated by physical considerations but 
most of them are phenomenological proposals for the desired energy density redshift |[l^-[|l4|. As a 
first step in the study of these models, the age of the Universe is calculated for several redshift laws 
of the energy density that resides in the dynamical scalar field [p^] , p^] . Observational consequences 



of an evolving A component decaying to matter and/or radiation have been studied in |17], [18|, 
obtaining severe constraints on such models. Another possibility, as the one that we consider here, 
is that the scalar field couples to matter only through gravitation. This kind of scalar field , with 
negative pressure and a time-varying, spatially fluctuating energy density, received the name of 



quintessence [19|. Its effects on the cosmic microwave background anisotropy are analysed in [19 
and [^] and phenomenological difficulties of quintessence have been studied in ||2l]] . Constraints on 
the equation of state of a quintessence-like component have been placed from observational data 



|22|. Recently, a potential for a cosmologically successful decaying A term has been constructed in 



m 



As we have discussed above, the behaviour of scalar fields is fundamental in understanding the 
evolution of the universe. In this paper we are interested in giving a general approach to the 
analysis of the cosmological evolution of scalar fields and to determine what kind of potentials 
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lead to a possible interpretation of the scalar field as quintessence and to a dominating energy 
density. However, we will not assume any kind of scale dependence for the potential nor impose 
any condition on which energy density dominates [|^, [p!^] , [p!^] , [p^]- p6[| . We will show that all 
model dependence is given only in terms of the quantity A = —V'lV^ where the prime denotes 
derivative with respect to the scalar field and its limiting behaviour at late times determines the 
evolution of the scalar field. 

Our starting point is a universe filled with a barotropic energy density, which can be either matter or 
radiation, and the energy density of a scalar field. The scalar field will have a self-interaction, given 
in terms of the scalar potential V{(\)\ but it will interact with all other fields only gravitationally. 
The barotropic fluid is described by an energy density and a pression with a standard equation 
of state = (7^ — l)/0^, where 7-^ = 1 for matter and 7^ = 4/3 for radiation. We do not make any 
hypothesis on which energy density dominates, that of the barotropic fluid or that of the scalar 
field. 

The equations to be solved, for a spatially flat Friedmann-Robertson-Walker (FRW) Universe, are 
then given by 

R = -^{p^+Pt + ^^) 

p = -2>H{p + p) (1) 

where H is the Hubble parameter, V{(p) is the scalar field potential, (/) = d(j)/dt, p (p) is the total 
energy density (pression) and we have taken 8ttG =1. It is useful to make a change of variables 
|2| X = <P/VQH, y = W/VSH and eqs.(|) b ecome 



XN 



VN 



-3x + ^Xy^ + ^x[2x^ + 77(1 -x^- y')] 

-^Xxy + ^y[2x' + 7,(1 - x' - y')] (2) 
Hn = -^H[j^{l-x^ -y^) + 2x^] 

where N is the logarithm of the scale factor a, N = ln{a), Jn = df/dN for / = x,y,H and 
A(A^) = —V'/V. Notice that all model dependence in eqs.@ is through the quantities A(A^) and 
the constant parameter 7^. Eqs.(|2|) must be supplemented by the Friedmann or constraint equation 
for a flat universe + + = 1 and they are valid for any scalar potential as long as the 
interaction between the scalar field and matter or radiation is gravitational only. This means that 
it is possible to separate the energy and pression densities into contributions from each component, 
i.e. p = p^ + P(j, and p = P7 + Pcf,, where p,f, {p^) is the energy density (pression) of the scalar 
field. We do not assume any equation of state for the scalar field. This is indeed necessary since 
one cannot fix the equation of state and the potential independently. For arbitrary potentials the 
equation of state for the scalar field p^ = (7,^ — l)/?,^ is determined once /?0,P0 have been obtained. 
Alternatively we can solve for x, y using eqs. (|2|) and the quantity 7,^ = (p^ +P<j>) / P4> = 2x^/ (x^ + y^) 
is, in general, time or scale dependent. 
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As a result of the dynamics, the scalar field will evolve to its minimum and if we do not wish to 
introduce any kind of unnatural constant or fine tuning problem, the minimum of the potential 
must have zero energy, i.e. V\min = V'\min = at (pmin- We will consider here only these kind of 
potentials. For finite (pmin the scalar field will naturally oscillate around its vacuum expectation 
value (v.e.v.). If the scalar field has a non zero mass or if the potential V admits a Taylor expansion 
around (pmin than, using the Hopital rule, one has limf^ool A| = oo and it will oscillate. On the other 
hand, if (pmin = oo then (p will not oscillate and A will approach either zero, a finite constant or 
infinity. The oscillating behaviour of i;^> or A is important in determining the cosmological evolution 
of x,y and = p^/rho = + ip' and we will show that any scalar field with a non- vanishing 
mass redshifts as matter field. 

Before solving eqs.(||) we define the useful cosmological acceleration parameter a and expansion 
rate parameter T. The acceleration parameter is defined as 

^ ^ p + 2,p ^ 37-2 

" ~ (37^ - 2)p 37^-2 ^ ' 

with 7 = (p + p)/ p. If Q = 1 then the acceleration of the universe is the same as that of the 
barotropic fluid and any deviation of a from one implies a different cosmological behaviour of 
the universe due to the contribution of the scalar field. A positive accelerating universe requires a 
negative a while for < a < 1 the acceleration of the universe is negative (deceleration) but smaller 
than that of the barotropic fluid. For a > 1 the deceleration is larger than for the barotropic fluid. 
In terms of the standard deceleration parameter q = one has a = ^^'^ ^ ™- terms of x,y 

one finds a = 1 — ^jh^iv^ ~ ^ ) — ^ ~ 3^1^ 3^_^22 • It is also useful to define the normalized 
equation of state parameter 

which gives the relative expansion rate of the universe with respect to the barotropic fluid. A F 
smaller than one means that the universe expands slower than the barotropic fluid and a F larger 
than one says that the universe expands faster due to the contribution of the scalar field. In our 
case a and F are not independent since F = 1 — (1 — q)(37^ — 2)/37^ = 1 — Vl^ "'''^ "''^ . 

A general analysis of eqs.(^) can be done by noting that, given the constant parameter 7^, all 
model dependence is through the quantity \{N). For an arbitrary potential V eqs.(|l]) or (Q) will 
be, in general, non-linear and there will be no analytic solutions. We can, of course, solve them 
numerically but we need to do it for each particular case and initial conditions separately. In order 
to have an understanding on the evolution of the scalar field we will study the asymptotic limit. It 
is useful to distinguish the different limiting cases for the cosmological relevant quantities x, y and 
= x^ + y"^ . Q.^ will either approach zero, one or a finite constant value. For $7^ the scalar 
field dilutes faster than ordinary matter or radiation and if fi^ — > 1 then the scalar dominates the 
energy density of the universe. When < cte < 1 then the scalar and barotropic energy 

density redshift at the same speed. Which behaviour will x,y,Q^ have, depends on A and on 7^. 
We will separate the analysis of eqs.(|2|) into three different behaviours of |A| at late times. In the 
first case we consider A a finite constant (or approaching one), A = c. Secondly we study the limit 
A — > 0. In the third case, we take A 00, which is the natural case if the v.e.v. of (p is finite but 
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we can have the same hmit for (p ^ oo. We divide in this case the analysis into an oscillating and 
a not oscillating |A| ^ oo. 

Eqs.(^) admit 5 different critical solutions for A constant |^6|. Here, we generalise these attractors 
for more complicated potentials that have a non-constant X{N). If x, y do not oscillate, since their 
value is constrained to |x| < 1, \y\ < 1, this implies that at late times they will approach a constant 
value, given by the attractor solutions of eqs.(||), and xn^Hn will vanish. Three of the 5 different 
critical solutions {xc = l,yc = 0), {xc = —l,yc = 0) and {xc = 0,yc = 0) are unstable (extreme) 
critical points. However, if A ^ oo then the critical point (xc = 0,yc = 0) becomes the asymptotic 
(stable) limit. The other two, depend on the value of X{N). 

For A^ > 37^ l2| one finds the "critical" values 
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for the quantities x,y, and an effective equation of state equivalent to that of the barotropic 
fluid (i.e. jtf, = {p(j, +P(f>)/P(j, = 77). In this limiting case the redshift of the barotropic fluid and 
the scalar field is the same. 

On the other hand if A^ < 6 then one obtains 

A^ 

= 1: 70 = y , (6) 

and the scalar energy density dominates the universe. If the scalar field has 7^ < 7^ then the 
solutions in eq.(^) are stable, the redshift of the scalar field is slower than that of the barotropic 
fluid. However, if = A^/3 > 7^, then the solution in eq.(^) is unstable and the scalar field ends 
up into the regime of the solution given in eq.(|5|). 

If A(A^) is constant then eqs.(||) or (^ are indeed solutions to xn = yN = Hn = but if A(A^) is 
not constant then the critical values in eqs.(^) and (^) solve xn = UN = Hn = only on single 
points, not an interval. This means that they are a solution to eqs.(§) only as an asymptotic limit 
and Xc{N),yc{N), X{N) are functions of A^. 

Let us now start with our first case, i.e. A constant. If A = —V'/V = c, the scalar potential has an 
exponential form, V = he~'^'^. This case has been extensively studied [18l,|^^-[26| and one finds 



critical (i.e. constant) points for x and y at late times. The value of x, y depends on the value of 
A = c and their solutions is given by eqs.(|5|) or (^). Since this case has been amply documented 



in the literature |18|,[25|-[^], we do not include its numerical analysis here. We show in fig.(2g) 
the behaviour of x and y for comparison with a double exponential potential {V = /le"'^'^'*) (see 
below). The cosmological parameters are a = r = lforA = c> ^37^, and a = ^7-7^, T = c^/37-y 
for A = c = -v/6- Note that in the first case is finite and even if it dominates the universe the 
acceleration and expansion of the universe is the same as for the barotropic fiuid. On the other 
hand, for A = c = \/6 one has Jl^ = 1 and the acceleration parameter a is in general different than 
one, it is negative if < 2 (assuming 7^ > 2/3 i.e. matter or radiation). In this case we could 
have interesting quintessence models. 
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For more complicated potentials, A is not a constant and its evolution determines that of x and 
y. The evolution of the scalar field leads to non-linear equations and critical points may exist but 
analytic solutions are either more difficult or impossible to obtain. However, the solutions given in 
eqs.(^) and (^) may give a good approximation on the limiting behaviour of x and y. 

Let us now consider the second case, i.e. A ^ 0. In this limit we can eliminate in eq.(§) the term 
proportional to A, and since —3 < H^/H < for all values of x,y and 7^, we have 

^ = -(3 + ^X0, ^ = -^>0. (7) 

X H y H 

From eq.(|^ we conclude that x will approach its minimum value (i.e. x ^ 0) while y will increase 
to its maximum value (i.e. y — > 1). In the asymptotic region with |x| <^ 1, |A| <C 1 one can solve 
eqs.(§) for x,y,H giving 



-3Af 
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with c, d integration constants. These solutions show that in the asymptotic region the scalar field 
dominates the energy density of the universe and the Hubble parameter goes to a constant value. 

In the limit A ^ 0, the first derivative of the potential approaches zero faster than the potential itself 
and examples of this kind of behaviour are given by potentials of the form V = V()(j)~'^,n > 0. The 
scalar field will dominate the energy of the universe leading to a " true" non- vanishing cosmological 
constant at late times with x — > 0, y — > 1, il,^ — > 1 and 7^ — > 0. The analytic solution to eq.(0) 
for x,y can be approximated by the expressions given in eqs.(^). However, x,y are no longer 
constant since Xc,yc depend on A(A^), which is itself not constant but these expressions are a good 
approximation at late times (see fig. Id) to the numerical results. 

The cosmological parameters are in this case a = — 3^^_2 and F = 0. This means that the 
acceleration of the universe is positive (since a < 0) and the expansion of the universe is exponential. 
In fig.(l) we show the behaviour of x,y,'y^ and for a potential V = Vb0~^. In fig.(la) we can 
see that the scalar field quickly dominates the universe behaviour and the Hubble parameter tends 
to a constant different than zero, i.e. the universe enters in an accelerated regime. This can also 
be seen in fig. (lb) where the acceleration parameter a is smaller than and 7,^, F are inflationary 
almost all the time. Fig.(lc) shows the behaviour of A(A^) for this case an din flg.(ld) we can see 
that the numerical solution has an asymptotic limit the solutions of eqs.(|^. 

As our flnal case we take the limit A — > 00, and we will separate this case into two different 
possibilities. The flrst one is when A approaches its limiting value without oscillating and the later 
case is when A does oscillate. 

In the non-oscillating case, in the region |A| ^ 1 the leading term of eqs.(|2|) is the one proportional 
to A if \y\, \x\ are not much smaller than one. In such a case the eqs. for XN-,yN are 

3 2 - /3 



XN = ]J-^y\, yN = -]J-^xyX. (9) 

The sign of xat is given by A and if it does not oscillate than x would reach its maximum value x = 1 
while y ^ for A > and x — > — 1, y — > for A < 0. However, in the region y ^ the other terms 
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Figure 1: Evolution of the universe filled with matter and a scalar field with V = Vo4>~^ ■ The initial conditions 
are xq ~ 0.1, yo — 0.5 and Ho — 1. In (la) we show (solid curve) quickly approaching 1; the Hubble parameter 
(dot -long dashed curve) tends to a finite constant, and for comparison, we have drawn Hm for a standard matter 
dominated universe (dotted line) ; notice that with this type of potential the difference in the rate of expansion with 
the standard model is remarkable. In (lb) we plot 70 (dot-long dashed curve), F (solid curve) and acceleration 
parameter a (dotted curve). In (Ic) A slowly evolves to zero with N. In (Id), the numerical x and y solutions are 
plotted (short and long dashed respectively) and compared to their attracting solutions eq. (^ , lower and upper solid 
lines respectively. 



6 



in eq.(§) become relevant and eq.(P) is no longer a good approximation. In the region |x| > |A| 
the evolution of x is given hy xn /x = —{?> + Hn / H) < and x will approach its minimum absolute 
value X — > like y. This region is the scaling region characterized by almost constant values of 
x,y and A. The end of the scaling region is when yN/u changes sign and becomes positive. This 
happens at ^3/2Ax + Hn/H ~ with Hjq/H ~ —87^/2. After a brief increase for y and x, they 
finally end up approaching the values given by the solution of eq.@ x ^ Xc,y ^ Vc and going to 
the extreme values oi x = y = 0. In this case Xc, yc are not really critical (constant) points since A is 
not constant. The kinetic and potential scalar energy will decrease faster than p, i.e. than ordinary 
matter or radiation. It is interesting to note that even though x, y approach zero, the equation of 
state of the scalar field becomes constant, i.e. 7^ = 2x^ j(x?' + y^) = cte because x^y decrease at 
the same velocity. This leads to 7^ approaching the value of 7^, i.e. the equation of state of the 
scalar field will be the same as that of the barotropic fluid (matter or radiation). Note in fig.(2h) 
that even though the equation of state of the scalar field approaches that of the barotropic fluid, 
its kinetic and potential energy decreases faster then that of the barotropic fluid, the reason being 
that 70 > 7^ at late times and the equality is only valid at i = 00. The cosmological parameters 
are a = T = 1 giving the same asymptotic behaviour for the universe with or without the scalar 
field. 

Examples of this kind of behaviour are given by potentials like V = e~°''^'^ , V = e~"'^'^ . In fig. (2) we 
show the behaviour of the dynamical variables and the cosmological parameters as a function of N 
for V = Ae~'^^ with A = —V'/V = ce'^. This potential gives the asymptotic limit for string moduli 
fields |2^. The solution of eqs.(^) shows that (f) ^ 00 minimises the potential and A ^ 00 at late 
times (fig. (2d)). The limiting values are x = y = = 0, as we can see in fig. (2a) and fig.(2h). 
In fig. (2a) we also show the evolution of the Hubble parameter in our model, as compared with 
the standard matter dominated case; we can see that the scalar field can influence the universe 
development only at early times. For small A^, the quantities x, y have a similar evolution as for a 
simple exponential potential. We exhibit these quantities in fig.(2g) for a typical case 7^ = 1 and 
A = 10. On the other hand, for N ^ 1, x and y are not constants and they approach the value given 
in eqs.(^), (A^ is larger than 87^), see fig.(2h). Fig. (2b) shows the behaviour of the 7^ parameter 
as X and y evolve and the effective total P parameter (c.f. eq.(Q)) for the "fiuid" composed by 
matter and the scalar field. Fig. (2c) is the acceleration parameter defined in eq.(^ Finally, fig.(2e) 
represents the phase space structure for (x, y) obtained with different initial conditions, where the 
final behaviour is amplified in fig.(2f). The plateau in the graph for A in fig. (2d), for approximately 
20 e-folds, corresponds to the scaling region, where x and y are constants (almost zero in this 
case), preceding the final evolution where the scalar field recovers a small quantity of kinetic and 
potential energy that finally go to zero. 

We conclude that such fields are not good candidates for quintessence (parameterizing a slow 
varying cosmological constant) and they do not play a significant role at late times unless they are 
produced at a late stage. We would like to emphasise that these results are completely general and 



leave out a great number of candidate fields such as string moduli |27|. The only condition we have 
used to derive these results is that A — > 00 without oscillating. 

We will now consider the case when |A| — > 00 but with an oscillating (p field. In this case the 
v.e.v. of (p is finite and without loss of generality we can take it to be zero. Around the minimum 
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Figure 2: Cosmological solution for V — A e~'^'^ , 7-y = 1 with initial conditions xq — 0.1, 3/0 = 0.5, Aq = 1 and 
Ho = 1. In (2a) the solid curve shows the evolution of O^, the dashed line is the numerical solution for H, while 
the dotted line is the comparison value of Hm for a standard matter dominated universe, as a function of A'^. In 
(2b) we show the evolution of 7^ (dot-long dashed curve) and of the effective total V parameter (solid curve). In 
(2c) the acceleration parameter a is displayed as a function of A^. In (2d) we have \{N). We show in (2e) the phase 
plane {x, y) for different initial conditions {xo, yo): (0.0, 1.0), (0.1, 0.9), (0.1, 0.7), (0.1, 0.5) and (0.1, 0.3); and the final 
evolution is amplified in (2f). (2g) compares the x and y evolution in this model (short and long-dashed curves 
respectively) with the solution obtained in the purely exponential potential V = he~'^'^ case (dotted and dot-dashed 
curves respectively), with A = c = 10, for small A'^. Finally, (2h) displays the asymptotic behaviour of x and y (curves 
are marked in the same way as in (2g)), showing the good agreement with the attracting solution eq.(^ (solid curve). 
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the potential can be expressed as a power series in and keeping only the leading term we have 
V = Vq (j)"" with n > and even since the potential must be bounded. The condition that 
V'\ 

min — requires n > 1 and a finite scalar mass requires n — 2. For this potential A — —n/(f) 
and it oscillates approaching a value |A| = oo, see fig. (3). As a first guess we could think that the 
limiting behaviour ol x,y is also given by eq.(^) and therefore tend to zero as A — > oo. However, 
this asymptotic behaviour is no longer a good approximation and we must solve the dynamical 
non- linear equations (see fig. (3) for the numerical solution of a quadratic potential). 

We will now determine under which conditions will either dominate (approach one), oscillate 
around a finite constant value or vanish. Since asymptotically H (x 1/t, a finite (7^ 1) requires 
that (j>, {(f))"'/'^ (X 1/t or equivalently that x,y are either constant or oscillate. We can thus write 
y = = B F"^^[G(t)], X = = A F2[G{t)] where Fi,F2 are arbitrary oscillating functions 
depending on a single argument G{t), and A,B are constants. The function G{t) is for the time 
being an unspecified function of t. Of course, -Fi and F2 are not independent since the functional 
dependence of (/> determines the functional dependence of 0, however this is not important at this 
stage. Taking the derivative w.r.t. N, we have y^ = ^2F\ ^nV and xn = G^x, where 
FiG = dFi/dG, = df/dN with f = x,y,G and i = 1,2. Since Fi,F2 are oscillating functions 
with a single argument G{t), we have that the average of < F^ >=< F^q > and the average of 
y^, is then 

2 

< >= Y < ^^^^ >' < '^^^^ > • 

In the asymptotic limit with oscillating and A — > 00 the evolution of x,y is given by eq.(^). 

Using eqs. 

(0) and d) we find that a potential V = Vq (j)^ may have finite values of x,y at late 

times and 

^ = ^. (") 

giving < 70 >= 2/(1 + -^f^) = 2n/(2 + n). Depending if 7^ is larger, equal or smaller than 7^, 
Qtf, will go to one, finite constant or zero respectively. Notice that this result is completely general 
and any massive scalar field redshifts at late times as matter fields with 7^ = 7^ = 1 (i.e. n = 2). 

In order to obtain the asymptotic solution analytically, we will solve the eqs.(||) in a region where 
only one component of the energy density dominates. This will be valid always in the asymptotic 
regime. However, we do not make any assumptions on which term dominates. If it is the barotropic 
fluid which dominates then the equation of state for p will have a parameter 7 = 7^, however if 
it is the scalar field that dominates than we take 7 =< 7</) > (average in time) since in this case 
'j^ oscillates. In this regime we can solve for the Hubble parameter and we get the standard form 
H(t) = 2/(37 1). To determine the evolution of the scalar field we have to solve its equation of 
motion (p + 3H(j) + ^^^^ = with V = Vq (p"" and n positive and even, i.e. 

4> + —<p + nl/o^"'^ = 0. (12) 
-ft 

For n = 2 the solution to eq. (^) is given in terms of Jm and the Bessel functions of the first and 
second kind respectively (j){z) = z~'"(2Vo)™'/^(ci Jm(-z) + C2Ym{z)) and 
= -z-™(2yo)('"-^)/2(clJ„^+l(z) + C2y„+i(z)) with ci,C2 Constants, m = (7 - 1/2), z = t^Wo 
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and we have used that d{z "^Km{z)) / dz = —z^Km+i with Km = JrmYm- Using these solutions 
we have 

y = rz'-^{kJUz)+Ym{z)) 

X = -rz^-"'{kJm+i{z)+Ym+i{z)) (13) 

with r = ^C2(2Vo)~, k = cijcx. A simple analytic expression can be obtained using the asymp- 
totic limit of the Bessel functions ^ ^2/-kz cos{z — ^(^™+^) ^ yj2/-KZ sin{z — Z!lH!B±li) 
for z 1 (i.e. t » 1). The amplitude of x and y in eq. (ll3|) in the limit z ^ \ goes as 
j; ~ y ~ Z"*^/^"™". A finite value of x,y requires m = 7 — 1/2 = 1/2 (i.e. 7 = 1). For 
7 > 1 (7 < 1) then X, y — > 1 (x, y — > 0) at large times. Furthermore, in the asymptotic limit 
7()i = 2[cos{z — ZLlHZB+li) _ k sin{z — ^(^'2+^) ))^/(l + A;^) is an oscillating function with an aver- 
age value < 7(^ >= 1. We can conclude therefore, that if the barotropic fluid has 7^ < 1, i.e. 
smaller than 7^, than = x^ -|- — > 0, but if 7^ > 1 =< 7^ > then the dominant energy 
density in the asymptotic regime will be that of the scalar field leading to Jl^ = 1. Finally, if 
7^ = 1 =< 7(^ > then the energy density of the scalar field dilutes as fast as the barotropic 
fluid and Vl^ tends to a constant finite value. The solutions in eq.(p!3|) for 7 = 1 can be given 
a completely analytic expression since in this case m = 1/2 and the Bessel functions take sim- 
ple form Jii2 = \f2pKZ sin(z), 1^1/2 = —\f2pKZ cos{z), J-^i2 = \plpKZ {sin{z)/z — cos{z)) and 
^3/2 = —\f2pKz {cos{z)/z + sin(z)). Putting these expressions into the definitions of x,y eq.([l^) 
we get 

y = yo sin{z) - {xq + —)cos{z), 

,sin(z) , , yo.,cos(z) , , 

X = -yo - cos{z)) + {xo + ^){^^ + sin{z)) (14) 

where the initial conditions are given by yo,X() at zq = 7r/2. In the limit z ^ 1 we have 0^ ~ 

The analytic solution in eq.(p^ agrees reasonably well with the one obtained by solving eqs.(0) 
numerically. This can be seen in fig. (3a), where we plot obtained numerically for 7^ = 1 in a 
dot-long dashed line and expression (pl|) in a solid line for a potential V = Vq (p"^. We also plot 
for 7^ = 4/3 and 7^ = 1/2 to illustrate the different asymptotic limits. For 7^ < 1 we have — > 0, 
7^ = 1 gives — > cte and for 7-^ = 4/3 we have 1. These different limits can be understood 

by noting that the average of < 7^ >= 1 and therefore, if 7^ > 7</)(77 < ^cj)) the barotropic fluid 
redshifts faster (slower) than the scalar field while for 7^ =< 7^ >= 1 both energy densities dilute 
at the same speed. The asymptotic value of in the limiting case A ^ 00 with an oscillating <p 
field depends on the value of 7^ and on the initial conditions xq = x(iVo),yo = y{NQ). Fig.(2b) 
shows the oscillating effective equations of state for the scalar field and the mixture matter/scalar 
field, and the resulting acceleration features of the universe a. 



We have obtained the analytic solution of eq.(|T^ for n = 2. This is clearly the simplest case since 
eq.(|l2|) is linear in cj) and it's derivatives. For n > 2 eq.(p^ becomes non-linear in cj) and no simple 
analytic solution exits. However, let us use in eq.(|l2|) the ansatz 



= (cicos(/3t2/'^) + C2sin{(3t^l'')) . (15) 
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It can be easily seen that this ansatz has the correct asymptotic behaviour (j)"'^'^ , 4>, H (x 1/t and 



oc 



if we wish to have x, y finite. The t exponents in eq.([l5|) are determined by solving eq.([l^) 
and this equation also imposes the conditions 7 = 2n/(2 + n) and /3 = jVon^C2~'^ {sin{l3t^/^) + 
k cos{l3t^/'^)Y-'^ . Notice that the value of 7 is precisely the value we obtained using general 
arguments only (c.f. eq.(ll)). Notice as well that the ansatz in eq.(^) is not a "complete" answer 



to eq.(|12|) since (3 is not a true constant. Only for n = 2, /? is indeed constant and the solution in 
dlD is the one we had previously obtained in terms of the Bessel functions (see eq.(|l3|)). For n 7^ 2 
we must take the average of f5 and work in the asymptotic region. Nevertheless, eq.(|l5|) gives a 
good analytic approximation to the numerical solution. In terms of eq.([T5|), x,y take the following 
expressions, 



k sin{at^l''') - cosiat^/'') + (sm(at2/") + k cos{at^''')) 



2/ti + k\ a 
y = yo (^sin{at^/'^) + k cos{at^/")y^'^ (16) 

with xq = ^1/3/2 P C2 7, yo = ^\/3Vo c^^"^ 7 and the initial conditions are taken at tg^" = it/ 2/3. 
Using eq.(|l6|) we obtain < y^ > j < x^ >= 2/n, at late times, as in eq.(^). We have, therefore, 
< 70 >= 2n/(2 + n) and 7 =< 7^ >, i.e. the ansatz in eq.(|l5|) is a "solution" to eq.(p^) only when 
the dominant energy density redshifts as fast as the scalar field. This is, of course, no surprise since 
we imposed on the ansatz, eq.(p^), the limit x,y — > cte. 

The cosmological parameters are in this case a = 1 — Sr^^ g^ ^2 F = 1 — '^^^ with 
70 = 2n/(2 + n). From these expression we conclude that when remains finite a = F = 1 since 
in this case 7^ = 7^ leading to the same behaviour of the universe with or without the contribution 
of the scalar field. However, when 17^ — > 1 then a = (870 — 2) 7(87^ ~ 2), F = 7^/7^ and one could 
have an accelerating universe if 7^ < 2 /8 which requieres that n < 1 and this is not acceptable since 
the first derivative of the scalar potential must vanish at the minimum. If Jl^ — > then clealrly 
a = F = 1 and the scalar field plays asymptotically no important role. 

To conclude this part of the analysis, we have established that if the initially dominant energy 
density component has a 7 parameter larger (smaller) than < 7^ >= 2n/(2 + n) then $7^ will 
approach one (zero). For n = 2 we have < 7^ >= 1, for n = 4 we have < 7^ >= 4/8. Since the 
condition V'\min = requires 1 < n we have that 7^ > 2/8 far all n and the scalar field will not 
give an accelerating universe. For n > 4 the energy density will decrease faster then radiation and 
since at late times the universe is matter dominated, only a scalar field with a non-vanishing mass 
could lead to a significant contribution to the energy density of the universe. However, since its 
redshifts goes as matter it is not a candidate for a cosmological constant but it could serve as dark 
matter. Fig.(8d) illustrates these characteristics of a power-law potential for the Hubble parameter 
H for V = Vo</'", n = 2,4, radiation and a matter dominated universe. 

To summarise and conclude, we have studied the cosmological evolution of the universe filled with a 
barotropic fluid and a scalar field with an arbitrary potential but only with gravitational interaction 
with all other fields. The analysis done is completely general and we do not assume any kind of scale 
or time dependence of the scalar potential nor any assumption on which energy density (barotropic 
or scalar) dominates. Our results are summarised in Table 1. We showed that all model dependence 
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Figure 3: Examples of evolution of the cosmological parameters for a universe filled with a perfect fluid (7^ — 
1,3/4, 1/2) and a scalar field with potential V = Vo^P . The initial conditions are xq = 0.1, j/o = 0.5, Ho = 1 and 
Vb = 37r'^/32. In fig. (3a) the numerical solution for Q.^ (dot-long dashed curve) with 7.^ = 1 is compared to the 
analytic expression Sl^ = + y'^ , calculated from ( p^ (solid curve). Also shown in the figure are Q,^ for 7-y =3/4 
(dot-short dashed curve) and 7-y = 1/2 (dashed curve). In (3b) the equation of state and acceleration parameters 
are displayed as a function of A*': 7^ (dot-long dashed line), V (solid line) and acceleration parameter (dotted line). 
In (3c) we plot A~^, in order to display the oscillating behaviour of \{N), as it approaches 00. In (3d) we plot the 
evolution of H for 2 different models. The solid curve represents the numerical solution for H, with V — Vo4>^ and 
7^ = 1 and is compared to H for a standard matter dominated universe (dotted curve). The dot -long dashed curve 
corresponds to our numerical solution for H, with V = Vo4>'^ and 7^, = 4/3, following a similar evolution as a standard 
radiation dominated model (dashed curve). 
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is given by A = —V'jV and 7^. Any scalar potential leads to one of the three different limiting cases 
of A: finite constant, zero or infinity. In the first case, 17^ approaches a finite constant (different 
than zero) depending on the value of A = c. For A — > we obtained x — > 0, y — > 1 with a constant 
Hubble parameter H and an accelerating universe. Finally, for A ^ 00 we concluded that if A does 
not oscillate x, y, fi^ — > and if A oscillates then all cases are possible (i.e. ^2,^ ^ 0, 1 or a finite 
constant) depending on the value of 7-^ and the power of the leading term in the scalar potential. 



m = -v'/v 


^,i> = P<i>IP 


70 


«((/)) 


r(<^) 




c =cte (> ^37^) 


'All 


77 


1 


1 




c =cte (< \/6) 


1 


3 


37-, -2 


A. 

377 




00 (no oscil.) 





77 


1 


1 




00 (oscil.) 




cte 
1 


^ 77) 
(= 77) 
^ « 77) 


1 
1 

37^-2 
37^,-2 


1 
1 

77 


Vb n > even 





1 





2 

377-2 





Vb > 



Table 1. In this table we show the asymptotic behaviour of 57^, 7^ the acceleration parameter a((^) — ^^~^2 ^^'^ 
the expansion rate parameter F = ^ for different limiting cases of In the last column we give an example of 

potential V{(j)) which satisfies this limit. 
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